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ABSTRACT 

We report our discovery of NLTT 20346 as an M5+M6 companion system to the tight binary 
(or triple) L dwarf 2MASS J0850359+105716. This nearby (~31 pc), widely separated (~7700 AU) 
quadruple system was identified through a cross-match of proper motion catalogs. Follow-up imaging 
and spectroscopy of NLTT 20346 revealed it to be a magnetically active M5+M6 binary with compo- 
nents separated by ^2" (50-80 AU). Optical spectroscopy of the components show only moderate Ha 
emission corresponding to a statistical age of ^5 - 7 Gyr for both M dwarfs. However NLTT 20346 
is associated with the XMM-Newton source J085018. 9+105644, and based on X-ray activity the age 
of NLTT 20346 is between 250-450 Myr. Strong Li absorption in the optical spectrum of 2MASS 
J0850+1057 indicates an upper age limit of 0.8 - 1.5 Gyr favoring the younger age for the primary. 
Using evolutionary models in combination with an adopted system age of 0.25-1.5 Gyr indicates a to- 
tal mass for 2MASS J0850+1057 of 0.07±0.02 M if it is a binary. NLTT 20346/2MASS J0850+1057 
joins a growing list of hierarchical systems containing brown dwarf binaries and is among the lowest 
binding energy associations found in the field. Formation simulations via gravitational fragmentation 
of massive extended disks have successfully produced a specific analog to this system. 
Subject headings: Astrometry- stars: low-mass- brown dwarfs- stars: fundamental parameters- 
binaries: wide- stars: individual (2MASS J08503593+1057156, NLTT 20346) 



1. INTRODUCTION 

Very low-mass stars and brown dwarfs (VLMs and 
BDs; M < 0.10 M ) are among the most populous con- 
stituents of the Galaxy yet their origins remain contro- 
versial. Whether they form in a manner similar to higher 
mass stars or require additional or com pletely different 
processes is currently under debate (see [B urgasser et al. 
[20071 [WTntworth et al.|2007| |Luhman et al.|2007| and ref- 
erences there-in). One important characteristic that in- 
dicates a difference in formation is their binary frequency 
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and multiplicity statistics (separation, mass ratio, eccen- 
tricities, etc). Contrary to higher mass stars, VLM bi- 
naries are found tightly bound (p < 20 AU), in small 
frequency (10-20%) and with high mass ratios (typically 
q~l - e.g. Bouy et al.|2003l |Close et al. " 



l2003l|Burgasser| 
il. 1120081). 



etal]|2003[ |Ahmic et al.||2007||Reid et al.||2008| ). 

Over the past decade, theoretical and observational 
studies have converged on several competing mechanisms 
to explain the formation of VLM stars and brown dwarfs. 
Among the most prominent are (1) magnetoturbulent 
or gravoturbulent fragmentation (Padoan fc Nordlund] 
20041 |Bate et al.|[2003[ |Ooodwin et al ||2004a|b[ |2006[ ); 
(2) ejection ot a pro tostclla r embryo from the na tal core 
jReipurth fc Clarke ||2001||Bate fc Bonnell||2005D; or (3) 



disk fragmentation ( Whitworth fc Stamatcllos j2"006| |Sta- 



matellos & Whitworth 2009a I. Characteristics ot multi 



2 



Faherty et al. 



pie systems predicted by each have been used to differ- 
entiate which mechanism is most probable for the pop- 
ulation. Early versions of the ejection model were fa- 
vored in part because they predicted that brown dwarf 
binaries should ha ve separations no greater than 10 AU 
This was supported through early 



2002) 



( Bate et al. 

observational studies. However recent findings of widely 
separated VLM systems (typically p > 100 AU) are con- 
tradictor y and shed doubt on the viability o f the ejection 



scenario (Luhman 2004: Allers et al. 2009 Close et al. 



20071 iLuhman et al ||2M9l JBilleres et al.JiWbSl |Caballero 
2007 |Artigau et al.|2007||Radigan et al.|2009| . Updated 
simulations have succ essfully created a, handful of wide 
brown dwarf binaries ( [Bate fc Bonnell||2005| |Bate||2009 1 
but the growing number in the field indicates a fraction 
too high to be explained by just this mechanism. The 
existence and frequency of wide substellar pairs may not 
discount one mechanism versus another; rather multiple 
formation mechanisms may be at work in the field. 

Recent work has shown that brown dwarfs widely sep- 
arated from nearby sta rs have a higher frequen c y of also 



being tight multiples (|Burgasser et al. 2005 Faherty 



et al.||2009 [). This h as been explored in the di s k Irag^ 
mentation models of Stamatellos & Whitworth (2009a I 
where similar high frequencies of brown dwarf multiples 
were produced when widely separated from a more mas- 
sive companion. In this article, we report the identifica- 
tion of another widely separated brown dwarf binary (or 
triple; see Burgasser et al 2010) companion, the NLTT 
20346 and 2MASS J0850+1057 comoving system. In Sec- 
tion 2, we review the discovery and observational data 
used to characterize the system. In Section 3, we discuss 
the astrometric details of the system including a con- 
taminating artifact that likely skewed previous parallax 
measurements for 2MASS J0850+1057. In Section 4, we 
analyze the details of the system including the activity 
and masses of each component. In Section 5, we dis- 
cuss its importance among the population of VLM wide 
companion systems. Results are summarized in Section 
6. 

2. OBSERVATIONS 
2.1. Targets 

2MASS J0850359+105716 (hereafter 2MASS 
J0850+1057) was first identified as a sin gle object 
in the Two Mi cron All Sky Survey (2M ASS; ]Skrutskie 
et al ||2006[ ) by |Kirkpatrick et al.| |l999[ ). It was subse- 
quently assigned to be the prototype for the L6 spectral 
su bclass before being res olved as an 7 16 bin ary system 



(20011. Bouy et al. (2008) confirmed 



by |Reid et aT] 

the binarity of 2MASS J0850+1057" through common 
proper motion detected in multiple epoch Hubble 
Space Telescope (HST) data. The magnitude difference 
between components of 2MASS J0850+1057 in HST 
images indicates that the secondary is a late-type L 



dwarf ( |Reid et al.||2001| |Bouy et al.||2008 ). Burgasser 
et al. 2010 report template fitting to the combined 
light spectrum of 2MASS J0850+1057 and conclude 
component spectral types of L7 and L6.5. They also 
speculate that the brighter but later-type primary may 
itself be a closely-separated, unresolved pair, making 
this system one of a hand ful of candidate brow n dwarf 



J0205 -1159, |Bouy et al.] [2005] Kelu 1, |Stumpf et al 
2009] , 

NLT T 20346 was i dentified as a high p roper motion 



star in 



Luyten| fll979|) and subsequently in Lepine et al 



( [2^02| . ^Euytenf ( |iy79 l assigned it a spectral class ol K 



using color information, but little else was known about 
NLTT 20346 prior to this work. 

2.2. ISPI Imaging 

We observed the field containing 2MASS J0850+1057 
and NLTT 20346 at the Cerro Tololo 4.0m Blanco tele- 
scope six times between 2008 March 09 and 2010 April 08 
as part of an ongoing astrometric program. The Infrared 
Side Port Imager (ISPI) was used and the target was 
obser ved through the J band filter (van der Bliek et al. 
2004 1 . ISPI has a nominal plate scale of 0.3" /pixel. Five 
images of 60 s exposures with 4 co-adds were obtained 
at each of the five dither positions. Observations were 
made under clear conditions with seeing ranging from 
-0.7-1.0". 

Dark frames and lights on/off dome flats were obtained 
at the start of each evening. The subsequent reduction 
procedures were based on the prescriptions put together 
by the ISPI tearrF] utilizing a combination of IRAF rou- 
tines as well a s publicly av ailable software packages: e.g. 
WCSTOOLS ( |Mhik][2002l ). J band flats were created by 
median combining the lights on and lights off images then 
subtracting the two. Bad pixel masks were created from 
a dome flat lights on and lights off image. Individual 
frames were flat-fielded and corrected for bad pixels with 
the resultant calibration images. All images were flipped 
to orient North up and East to the left using the IRAF 
routine osiris in the cirred package. Finally, the IRAF 
routine xdimsum was used to perform sky subtractions 
and mask resultant holes from bright stars. 



NLTT ^346B 



£ " Sep 



NLJT20346A 



Fig. 1. — SpeX J band image of the 6x6 arcsccond field around 
NLTT 20346 where we resolved the system into two M dwarf com- 
ponents. The angular separation is ~ 2 " (or 50-80 AU) and the 
J band magnitude difference is 1.13±0.02. 

In close examination of an ISPI image from 2009 
November 30 (seeing conditions ~0.7"), we noted an 
elongated point spread function (PSF) for NLTT 20346 



triples (e.g. Gliese 569, Simon et al 



2006] DENIS 



http : I /www. ctio.noao.edu/instruments/ir_instruments/ispi/ 
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in J band after a Is exposure indicating a visual binary. 
Subsequent images taken with the Cont — 203 filter (and 
J, H, K filters on the SpeX guider camera-see section 
2.3.2 belowWesolved a companion separated by ~2.0 " 
(see Figure IT]) . 
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Fig. 2. — The optical spectrum using MagE of the primary (top 
plot) and secondary (bottom plot) components of the NLTT 20346 
system. Over-plotted is the template for an active M5 (top plot) 
and M6 (bottom plot) from ^Bochanski et al. (20071 normalized at 
7400 A (dotted line). The inset shows Ha (6563 A) emission and 
a lack of Li (6708 A) absorption in both components. 



2.3. Spectroscopy 
2.3.1. MagE 

Optical spectra for the components of the NLTT 20346 
visual binary were obt ained with the Mage llan Echellette 
Spectrograph (MagE; |Marshall et al.||2008| on the 6.5m 
Clay Telescope at Las Campanas Observatory on 2008 
November 26 and 2010 March 9 (see Figure M. MagE 
is a cross-dispersed optical spectrograph, covering 3,000 
to 10,000 A at medium resolution (R ~ 4,100). Our 
observations employed a 0.7" slit aligned at the parallac- 
tic angle. Observations were made under clear conditions 
with an average seeing of ~0.7". The separation between 
the primary and secondary is ^2" therefore minimal con- 
tamination (< 1%) was expected. A 300s integration was 
obtained for each component followed by a ThAr lamp 




1.5 2.0 
Wavelength (nm) 



I 

ft 


' 1* H,0 








\ 






V 

J 

1 




FeH 

1 H,0 










CO 


1 




1 I 1 





1.6 2.0 
Wavelength (nm) 



Fig. 3. — The near-IR spectrum using SpeX of the primary (top 
plot) and secondary (bottom plot) components of the NLTT 20346 
system. Over-plotted as a dotted line is the M4 optical standard Gl 
213 on the primary component (top) and the M5 optical standard 
Gl 51 on the secondary component (bottom). 

spectrum for wavelength calibration. The spectrophoto- 
metric standard GD 108 or EG 274 was observed for flux 
calibration. Ten Xe-flash and Quartz lamp flats as well 
as twilight flats were taken at the start of the evening for 
pixel response calibration. The data were redu ced using 
the MagE Spe ctral Extractor pipeline (MASE; Bochan- 
ski et al.]|2009 ) which incorporates flat fielding, sky sub- 
traction and flux calibration IDL routines. 
2.3.2. Near-Infrared Spectroscopy and Imaging with SPEX 

We obtained low-resolution near-infrared (NIR) spec- 
troscopy for both the primary and secondary com ponents 
of NLTT 20346 using the SpeX spectrograph (Rayner 



et al. 2003 ) mounted on the 3m NASA Infrared Telesco 
Facility (1RTF) on 2009 December 3 UT (see Figure [3} . 
We used the spectrograph in prism mode with the 0'.'5 slit 
aligned to the parallactic angle. This resulted in R = A 
/ AA « 120 spectral data over the wavelength range in 
0.7-2.5 /im. Conditions were clear and the seeing was 
0746 at K . The stars were separated sufficiently that we 
were able to obtain 12 individual exposure times of 30 
seconds (total exposure time of 6 minutes) in an ABBA 
dither pattern along the slit. Immediately after the sci- 
ence observations we observed the AOV star HD 43607 at 
a similar airmass for telluric corrections and flux calibra- 
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tion. Internal flat-field and Ar arc lamp exposures were 
acquired for pixel response and wavelength calibration, 
respectively. All dat a were reduced using the SpeXtool 



packa ge version 3.4 (Cushing et 
2003 1 using standard settings. 



al. 2004 Vacca et al. 



In addition to NIR spectroscopy, we used the SpeX 
guider camera to image NLTT 20346 and resolve the two 
components with a S/N of > 50 for each. The resolution 
is Of! 12 per pixel with a 60" x 60" field of view. We used 
the JHK filters and the camera oriented with North up 
and East to the left. Exposure times were 5 seconds each 
in an ABBA dither pattern, which was then shifted 10" 
to the east and repeated (see Figure [3]). A sky image was 
created from the 8 images then subtracted from individ- 
ual frames before relative photometry was performed. 
3. ASTROMETRIC DETAILS 

3.1. Contaminant of 2MASS J0850+1057 

2MASS J 0850+1057 was tar geted in the astrometric 



programs of |Vrba et al] ( |2004[ ) and |Dahn et al | ( |2002 1 
resulting i n two published b ut discrepant parallax val- 
ues. The Vrba et al. (2004) work used a near-IR im- 
ager with either J or H filters (H in th e case of 2MASS 
J0850+1057) and the |Dahn et~aL] ( |2002j) work used an op- 
tical CCD with a wide — I interference filter for astromet- 
ric measurements. There were seven objects in common 
between the two programs, five of which had parallax 
measurements that matched better than 4 mas. 2MASS 
J0850+1057 was discrepant by 12.9±5.5 mas (CCD-PR) 
resulting in a 2.4c difference in parallactic angle and a 
3.4(7 differe nce in proper motion position angle. |Vrba| 
et al. (2004) investigated whether orbital motion between 
the resolved components of the binary could account for 
the disagreement but ruled it out due to the short time- 
span between the two programs compared to the pre- 
dicted ~51 year orbit. 

We noticed on recent ISPI images that there was a sec- 
ond fainter object resolved in the J band ~ 4.0" from the 
L dwarf that could have ske wed prior astrome tric mea- 
surements (see Figure H). The | Vrba et al.|||2004fr parallax 



measurement was based on 1.86 year s of data, 13 i mages 



(2002) 



and a mean epoch of 2001.791. The Dahn et al 
parallax measurement was based on 3.3 years of data 
and 30 images. Since no mean epoch is given, we assume 
observations were conducted betwee n 1999-2002 s hortly 
after the discovery and before the Dahn et al. (2002) 
publication. Both astrometric programs have compara- 
ble resolutions to ISPI. Tracing the position of 2 MASS 
J0850+1057 back to the |Vrba et all j2004| and [Dahn 
et al. ( 2002 ) mean positions from the ISPI position using 
our updated proper motion vector revealed that the bi- 
nary L dwarf would have been blended with the contami- 
nant at the time of the cited astrometric programs. Over 
the 1.86-3.3 years of the programs, the significant proper 
motion of 2MASS J0850+1057 would have changed the 
centroid shape as its separation with the contaminant 
increased. This potential "elongation" of the centroid 
would have skewed the position of 2MASS J0850+1057 
in the frames used to calculate the parallax and proper 
motion. 

A digital sky survey image from 1995 Nove mber 17 in 



J ban d as well as a Sloan Digital Sky Survey (York et al 
2000 1 image from 2005 November 06 in the sloan i band 
also show the point source. The second object appears 
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Fig. 4. — ISPI J band image of the 10x10 arcsecond field around 
the target 2MASS J0850+1057. The point source at the far left 
is ~ 2 magnitudes fainter at J but clearly resolved at the 2009 
ISPI position. Over-pl otted with filled cir cles are the mean epoch 
positions from 2MASS, |Dahn et al 1(2002), and | Vrba et al,|j2004). 
Tracing back the ISPI position to tne JVrba et al.| ( |20U4p ana fL>ahn| 
|et al.| ( |2002[ ) positions using the proper motion calculated in tnis 
work reveals that 2MASS J0850+1057 would have been blended 
with the contaminant. Over the lifetime of the two astrometric 
programs, the significant proper motion would lead to an "elon- 
gation" of the PSF which is a likely explanation for discrepant 
parallax measurements between the two programs. 

to be 2 magnitudes fainter in J but shows no appreciable 
motion and we conclude that it is unrelated to the brown 
dwarf. We propose that this contaminating source is the 
most probable cause for the discrepant parallax results. 

3.2. Proper Motion 

2MASS J0850+1057 was imaged as part of the Brown 
Dwarf Kinematics Project (BDKP) parallax program 
which targets nearby brown dwarfs. We used the 
Carnegie Astrometric Planet Search software (from here- 
on ATPa) to extract all point sources in the 6 epochs 
of IS PI data and solve for the parallax and proper mo- 
tion (Boss et al. 2009). The highest quality image was 
used as the template which all other images were trans- 
formed. In the 2MASS J0850+1057 field, there were 
25 well-behaved (elongated, saturated, spurious sources 
were removed) reference stars between all epochs. Us- 
ing these sources, a linear transformation was applied to 
each epoch point source catalog to constrain the field ro- 
tation, plate scale and match all reference sources to the 
template. Higher order transformations were tested but 
demonstrated negligible difference from linear solutions. 
The apparent trajectory of each star was fit to a standard 
astrometric model included in the ATPa software. The 
algorithms follow the astro metric solution prescrip tions 



laid out in the Hi pparcos ([Perryrnan et 



Tycho Catalogue (|Hog et al.||2000|) description 
corrected the data 



from apparent parallax to absolute 




5 http : //www.rssd.esa.int/S A/HIPP ARCOS/docs/vollJl.pdf 
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TABLE 1 
Details on Astrometry 



Name 
(1) 


RA 

IZjpOCfl zuuu 
(2) 


DEC 
H/pocn zuuu 
(3) 


fJ>a 
ii —l 

y r 

(4) 


W 
« w — l 

(5) 


^rad 

km s 
(6) 


Distance 
pc 
(7) 


Reference 
(8) 


NLTT 20346 a 
NLTT 20346 
NLTT 20346 


08 50 19.18 


+10 56 43.69 


-0.172 ± 0.008 b 
-0.179 ± 0.010 
-0.177 ± 0.019 


-0.050 ± 0.006 b 
-0.052 ± 0.010 
-0.035 ± 0.019 


26 ± 9 C 


31 ± 7 d 


1 

2 
3 


NLTT 20346A 
NLTT 20346B 






-0.178 ± 0.033 c 
-0.166 ± 0.033° 


-0.050 ± 0.032 c 
-0.050 ± 0.032 c 


25 ± 7 
26+7 


28 + 5 
33 + 5 


1 
1 


2MASS J0850+1057 
2MASS J0850+1057 
2MASS J0850+1057 


08 50 35.93 


+10 57 15.62 


-0.144 ± 0.006 
-0.143 ± 0.006 
-0.142 ± 0.002 


-0.038 ± 0.006 
-0.020 ± 0.003 
-0.008 ± 0.002 




29 + 7 
38 + 6 
25.6 ± 2.3 


1 

1 
5 





References. — Refs-l=This paper 2= jLepine et al.||2002) 3= |Luyten|jl979( 4= |Vrba et al.| ( |2004^ 5= |Dahn et al.H2002) 
a NLTT 20346 consists of two M dwarf components separated by ~2" or 50 - 80 AU. 
b Based on the catalog positions of the blended combined light source over a 59 year baseline 
c Based on the mean RV value of NLTT 20346 A and B. 

Based on the mean value from the spectrophotomctric distance to the components. 
c Based on the positions from resolved SDSS and ISPI images with a ~4 year baseline. 



parall a x follo wing the same procedure described in |Vrba| 



et al. (2004). A full detailed description of the astro 
metric pipeline used to solve for the parallax and proper 
motion of 2MASS J0850+1057 is provided in Fahcrty et 
al. (in prep). As noted in Vrba et al (2004) the proper 
motion and parallax of the system will not be effected 
by the orbital motion since the images were taken over 
a ~ 2 year period whereas the orbit is ~ 50 years with 
components separated by 0.16". 

To determine the new proper motion for NLTT 20346 
we used positions available through USNO-A2.0, USNO- 
B1.0, GSC 2.2, 2MASS, and an ISPI image taken on 08 
April 2010 ([Monet et al. | |2003[ |Monet||1998|Lasker et ah 
20081 ICutriet al. 2003). Thetotal baseline between the 



first and final epoch used was ~59 years. NLTT 20346 
is not resolved in any of the epochs as the 2" separation, 
seeing conditions and plate scale of each detector caused 
the A and B components to appear blended. We simu- 
lated the PSF for the A and B component of NLTT 20346 
on each of the detectors to estimate the uncertainty due 
to blending. We added this additional uncertainty in 
quadrature to the catalog uncertainties and used these 
in the proper motion measurement. The absolute astro- 
metric position from each catalog (with corresponding 
epochs) were used to solve for the proper motion using 
a least squares weighted solution. There is an SDSS im- 
age taken on 06 January 2006 and an ISPI image taken 
on 27 February 2010 with sub-arcsecond seeing where 
both components are resolved. Using these two epochs 
with a ^4 year baseline, we calculated the proper motion 
of each component. These resolved values are consistent 
with the motion calculated from the blended source along 
a 59 year baseline. This indicates that orbital motion is 
not effecting the total motion of the system. 

For NLTT 20346 we calculated fi a =-172±8 mas yr" 1 
and /i,5=-50±6 mas yr" 1 (using all catalog positions) and 
for 2MASS J0850+1057 we calculated /^ Q =-144±6 mas 
yr" 1 and ^ 5 =-38±6 mas yr" 1 . NLTT 20346 has two 
published prope r motion value s in the LSPM-N and New 

catalogues and, as stated 
+1057 has prev ious proper 



1979) 



Luyten (NLTT- |Luyten 
in section 3.1, 2MASS . 

motion valu es reported in |Vrba et al. (20041 and Dahn 



consistent within 2a of previous published results and 
using the new values , both the fi a and jj,g values for 
the potential companions are within 2a of each other. 
Along with the predicted distance measurement of the 
primary (see section 4.1 below), and the new parallax 
measurement for 2MASS J0850+1057 of 35±8 mas or 
29±7 pc, this system is a strong wide companion common 
proper motion candidate. 

3.3. Likelihood of Companionship 

NLTT 20346 was identified as a potential compan- 
ion to 2MASS J0850+1057 through a common proper 
motion search of the Brown Dwarf K inematics Project 
(BDKP) catalog ( |Faherty et~aL]|2009[ ) and the Lepine- 
Shara Pr oper Motion INorth (LSPM-N) and Hip parcos 
catalogs ( Lepine et al.|200"2 Perryman et al.|1997 ). Sev- 
eral other syste ms were detected and detailed analysis 
was presented in Faherty et al. (2009). In that work, an 



angular separation ot up to 10 aremmutes and a proper 
motion match criterion of better than 2a in both right as- 
cension (RA) and declination (DEC) between the system 
components was required to determine common proper 
motion candidates. The average uncertainty for objects 
in the BDKP catalog is 15 mas yr -1 (in both directions) 
so proper motion agreement was typically required to be 
< 30 mas yr -1 bet ween the stellar comp anion and ultra- 

(2010J) also required a 
of better than 2er 



cool dwarf (UCD). Faherty et al 



et al. (2002) (see Table [I]). Our astrometric results are 



distance match between components 
or typically better than 10 pc. 

The sy stem containing NLTT 20346 was not investi- 
gated in |Faherty et al" (20101 because proper motion 
components were slightly outside of the 2a requirement. 
However, after follow-up imaging and re-analysis of the 
astrometry of both components we found strong evidence 
for companionship (see Section 3.2). 

To quantify the probability that NLTT 20346 might be 
a chance alignment with 2MASS J0850+1057, we ran a 
Monte Carlo simulation of all stars in the LSPM-N and 
Hipparcos catalogs that shared a common proper motion, 
but not necessarily dist ance or position, with the brown 
dwarf (to within 2<r-see Faherty et al. 2009 for details). 
There were 156 stars in the Hipparcos catalog and 632 
stars in the LSPM-N catalog with matching proper mo- 
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tion components. After 10000 iterations we found the 
likelihood that NLTT 20346 is a chance coincidence with 
2MASS J0850+1057 (at an angular seperation of 248 ") 
is < 1.0%. 

4. SYSTEM PARAMETERS 

4.1. Binary Components of NLTT 20346 

From the MagE data for NLTT 20346 we classified the 
brighter component as an M5±0.5 and the fainter as an 
M6±0.5 by c omparing to active M dw arf templates (see 



Bochanski et al. (2007). SpeX imaging of 



Figure k2[) in 

NLTT 20346 allowed us to measure the relative mag- 
nitude difference between components (A J=1.13±0.02) 
and subsequent J band magnitudes of 11.61± 0.04 and 
12.74± 0.08 for the M5 and M6 respectively. From the 
SpeX spectrum, we classify the brighter of the two com- 
ponents as an M4±l and the fainter as an M5±l by com- 
paring to the near-IR spectra of the M4 optical standard 
Gl 213 and the M5 optical standard Gl 51 respectively 
(see Figure [3]). This is consistent within uncertainties 
with the optical spectral types. In the absence of a par- 
allax measurement we computed a spectrophotometric 
distance to each component using the spectral type, indi- 
vidual J band magnitudes, and the spectrophotometric 
J band relation from Golimowski et al (2010 in prep). 
We calculated component distances of 28±5 pc and 33 
± 5 for the M5 and M6 respectively resulting in a mean 
distance value to NLTT 20346 of 31±7 pc. 

4.2. Activity and Age of NLTT 20346 

The optical spectra of both components of NLTT 
20346 show moderate Ha emission. We measured an Ha 
equivalent width of 4.20±0.06 A and 3.64±0.08 A from 
the MagE data for the M5 and M6, respectively. The 
upper Balmer series showing H/3 through H<5 as well as 
Ca II K and H+He are also seen in emission in both 
components (see Figure [2b. Combinin g the Ha equiva- 
lent w idth with the \ parameter from Walkowicz et al. 
(2004) yields the log(L# a /L& i), a me tric of magnetic ac- 



tivity and a statistical proxy for age. |West et al] (2008[ ) 



have found mean log(Ljj a /L(, ;) values tor active Ma and 
M6 dwarfs within 100 pc of the Sun (see Figure [5]) of - 
3.9±0.2 and -4.0±0.3, respectively. NLTT 20346A and B 
have log(L ffQ /L bo; ) values of -3.94±0.02 and -4.20±0.02 
respectively , both within la of typical active M dwarfs. 

12008| estimate an activity lifetime for M5 
of 7. 0±0.5 Gyr . Using the age activity 



West et al. 
and Mb dwarfs 



relation from [West et al. ( 2009jn yields consistent com- 
ponent ages of 6.3±1.0 and 6.5:£1.0 Gyr for the M5 and 
M6 respectively. 

The NLTT 20346 system is an X-ray source and the 
X-ray activity indicates a younger age than that calcu- 
lated from the Ha activity. 2XMMiJ085018.9+105644 
(from the X-ray Multi Mirror Mission(XMM)-Newton 
2nd Incremental Source Catalogue) is 2" from the po- 
sition of the source. We computed the X-ray flux (fx) 
by combining the 0.2-0.5 keV, 0.5-1.0 keV, and 1.0-2.0 
keV ranges in order to compare to equivalent X-ray de- 
tections of M dwar fs in the Rontgen Satellite (ROSAT; 



Voges et al. 1999|) catalog (~ 0.2-2.4 keV ) We used 



the i x /ij relation defined in Agiieros et al. ( |2009 ): log 



6 Equation 3.1 in |West et al. ] |2009| | contains an error. The 
functional form should be log(Ljf a /Lj, ;)= pESm ~ b 
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Fig. 5. — Distribution of loafL^^/L i,,,)) for active M5 ( top plot) 
and M6 dwarfs (bottom plot) in the IWest et al.l (|200Sb spectro- 
scopic sample within 100 pc of the Sun. 1NL,1"1 ^U34bA (top) and 
NLTT 20346B (bottom) have values of -3.94±0.02 and -4.20±0.02 
respectively (marked by a five point star), within lcr of typical 
active objects at similar temperatures. 

(f a /f^=log(/ x )-r-0.4.J+6.30 to compute log for 
NLTT 20346 of -2.5. Note that the 2MASS J band mag- 
nitude for the combined light of the system was used to 
calculate log (f x /fj) since the components were not re- 
solved in the XMM-Newton data. Figure [6] compares 
this estimate to Hyades, Pleiades, Young (< 300 Myr) 
and thin-disk field M stars flShkolnik et al.|2009|) . The ob- 

! IT 



jects labeled Field Ms are from the Himsch et al.| ( 1999| ) 
study. Several of these objects de monstrate large X-ray 
flux although they are regarded in |Shkolnik et al. (20091 
as normal field objects (typically > 1 Gyr). P ossible ex- 
planations include flaring M dwarfs (Hilton et al. 2010 
calculate a flare duty cycle for late M dwarfs of as much 
as 3%), or undetected spectroscopic binaries (Shkolnik 
let al.l [20091 found a 16% contamination rate of SBs in 
t heir sample) 



Preibisch & Feigelson ( 2005 1 investigated the age- 
activity trend for low mass stars. Comparing X-ray ac- 
tivity for stars in the mass range 0.1 <M Q <0.4 belonging 
to nearby associations or clusters with ages spanning 1- 
100 Myr (comparison samples were ONC, NGC 2264, 
Chamaeleon, Pleiades, and Hyades) they found only a 
mild decrease in X-ray luminosity over the 1-100 Myr 
range. However for the same mass range they found 
that significant decay was detected over longer timescales 
(such as from Hyades age objects to Field-several Gyr- 
age objects). Unfortunately no relation is derived for 
obtaining an age for the slightly older population of M 
dwarfs based on X-ray luminosity. On Figure [6[ the X- 
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Fig. 6.— Comparison of the X-ray activity of NLTT 20346 to 
M dwarfs i n the Pleiades, Hyade s, young and older field M dwarfs 
(data from Shkolnik ct al. 2009). NLTT 20346 is marked by the 
five point star. JNote tnat the log (fz/fj) for NLTT 20346 is based 
on the combined light in J band and X-ray as we do not have 
resolved X-ray data for the components. 

ray activity of NLTT 20346 places it with the most active 
Hyades stars or the least active Pleiades stars showing 
similar colors. It also has X-ray activit y levels compa- 
rable t o the young M dwarf sample from Shkolnik et al. 
( 2009 ) where the estimated age of the sample was < 300 
Myr (with objects well above the X-ray activity level 
of the Hyades having ages < 150 Myr). To estimate the 
age from the X-ray activity level we assume a Skumanich 
like decay from the mean log (f x /ij) of the Pleiades (-2.3 
±0.3) with an assumed age of ~ 120 Myr. The resulting 
age range for NLTT 20346 is 250-450 Myr consistent with 
the young field M dwarf sample for which it has compa- 
rable activity le vels. We note t hat th ere are a handful of 
M dwarfs in the Hiinsch et al. ( 1999 ) sample which show 
comparable activity levels indicating that we cannot rule 
out that NLTT 20346 is a flare star or spectroscopic bi- 
nary (explanations that could account for the high X-ray 
levels) . 

The ^2" or 50-80 AU separation of the components 
of NLTT 20346 is safely outside the distance where tidal 
interactions would affect stellar rotation and hence the 



chromosp heric age of the system ( Meibom et al. 2007). 



However, |Armitage fc Clarke ( 1996 1 find that binary sys 
tems with separations <100 AU can suffer accelerated 
mass-accretion from star-disk interactions early in their 
lifetimes which could lead to variations in the rotation 
rate. The NLTT 20346 system is close to the border 
where disk interactions become significant, so the X-ray 
age above might be skewed to a slightly younger age. 

We cross-correlated the MagE spectra of NLTT 
20346A and B with the M5 and M6 templates from 



IBochanski et al. ( 2007 ) using the xcorl pac kage in IDL 
QMohanty fc Basri|2003 [ |West fc Basri||2009[) to measure 
consistent radial velocities of 25±7 km/s and 26±7 km/s 



respectively. Combining the mean RV value of 26±9 
km/s with the proper motion, position, and spectro- 
photometric distance leads to (U,V,W) values of (21±5,- 
10±4,-1±4) km/s (corrected for solar motion with U 
positive toward the Galactic center). While, individual 
space motions can not be used to date objects, the calcu- 
lated Galactic space motion for NLTT 20346 is consistent 



with thin disk membership (age < ^ 3 Gyr), in agreement 
with the chr omospheric diagnostics ( Eggen fc Iben|1989~ 
Eggenfl989 l 



Discrepancies between various age diagnostics are not 
surprising. Previous studies have shown similar in- 
consistencies when determining ages for individual ob- 
jects using different techniques (e.g. isochronal ages, 
gyrochronol gy, chromospehri c activity, Lithium abun- 



dances; see Soderblom 2010b. Some possible explana- 



tions for the differing age values of this system are that 
(1) we could have observed the components of NLTT 
20346 when the system was in a low state as variations 
of more than 30% have been seen in Ha measurements 
of M dwarfs (e.g. |Berger et aT1|2009[). (2) Preibisch & 



Feigelson (2005) show the X-ray levels of objects < 500 
Myr versus field age objects with median ages of a few 
Gyr; however very little is known of what to expect for 
intermediary or juvenile age objects. In this work we 
favor the younger age for the system as this shows con- 
sistency with the presence of strong Lith ium absorption 
in 2MASS J0850+1057(see section [o]), the X-ray ac- 
tivity in NLTT 20346, and the kinematics of the system. 

4.3. Age of 2MASS J0850+1057 
The optic a l spec trum of 2MASS J0850+1057 from 



Reid et al. 



( |2001[ ) exhibits deep Li absorption (15.2 
A equivalent width) indicating component masses below 
~0.06 and an upper age limit betwee n 0.80 - 1.5 Gyr 
( |Rebolo et al.|1992||Magazzu et al.|1993[ ). The strong Li 
absorption is consistent with the X-ray activity and sub- 
sequent younger age of NLTT 20346. Coupled with the 
matching proper motion components and distance mea- 
surements further confirms the likelihood of companion- 
ship. 

We choose two age ranges for further analysis of this 
system: (1) 0.25 - 1.5 Gyr based on a lower limit con- 
strained by the X-ray activity of the primary and an up- 
per limit based on the Li detection in the secondary and 
(2) 5 - 7 Gyr based on the Ha activity of the primary. 
4.4. Component Masses 

We determined masses for all components of the new 
system using the two age ranges discussed above (0.25- 
1.5 Gyr and 5 - 7 Gyr), Lb i values (described in this 



( 1997 ) evolutionary mod- 



values for 2MASS 



section) , and th ejBurrows et al. 
els (see Figure [7| . To comput 
J0850+1057 we combined our new parallax measure- 
ment with the resolved near-IR photometry and spec- 
tral type components reported in Burgasser et al (2010), 
and the bolo metric correction calculated from |Dupuy| 
& Liu (2009). We determine log(Lfc //L su „) values of 



4.43±0.20 and -4.82±0.20 for the primary and secondary 
respectively. Using evolutionary models and the younger 
age range of 0.25-1.5 Gyr, the calculated component 
masses for 2MASS J0850+1057 are 0.04±0.02 M@ and 
0.035±0.015 M yielding a total brown dwarf binary 
mass (M tot ) of 0.075± 0.025 M Q . Using the older age of 
5-7 Gyr the component masses are significantly higher; 
0.076±0.02 M and 0.072±0.03 M Q for the primary and 
secondary respectively yielding a total brown dwarf bi- 
nary mass (M tot ) of 0.148± 0.040 M . Table [2] lists all 
s ystem parameters. 



Konopacky et al. (2010) calculated the total mass of 



this system to be 0.2± 0.2 M Q based on Keck AO obser- 
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TABLE 2 
Details on System Components 



Parameter 2MASS J0850+1057 NLTT 20346 Reference 



Combined Spectral Type 


L6 


M5 


1,2,5 


Est. Component Types 


L7+L6.5 


M5+M6 


3,5 


2MASS J 


16.465 ± 0.113 


11.282±0.023 


4 


2MASS K 


14.473± 0.066 


10.407±0.020 


4 


A Jmko 




1.13±0.02 


5 


A K MKO 


0.95±0.02 




3 


Primary log(L(, oi /L ) a 


-4.43±0.20 


-2.46±0.09 


5.7 


Secondary log(Lj, oi /L0) 
Primary Mass (M Q ) a 


-4.82±0.20 


-2.69±0.09 


5.7 


0.04±0.02 


0.16 ±0.01 


5 


Secondary Mass (Mq) 


0.035±0.015 


0.13±0.01 


5 


System Age (Gyr) 


0.25-1.5 


0.25-1.5 


5 



References. — Refs-1 = |Kirkpatrick et al.| jl999|; 
< |2000|> ;3=Burgasser et al (2 010); 4= |CJutri et al.| $2 
6= |K.onopacky et al.| ( |2010t ; 7= |Reid fc Hawley| l |2005| l 
a Assuming this component is a single; sec Burgasscr ct al. (2010). 



Kirkpatrick et al.| 



2003 ~ his paper; 



-2.0 ~ 




-0.5 0.0 
log(age)[Gyr] 



Fig. 7. — A plot of the |Burrows et aT7| ( |1997[ l evolutionary models 
with parameters (age and 1jj, ; ) lor both the primary and secondary 
of each component of the quadruple system indicated with labeled 
boxes. Masses are shown derived for the favored age range of 0.25 
- 1.5 Gyr as well as the He* age range of 5 - 7 Gyr. 

vations and orbital fitting. The large uncertainty of the 
dynamical mass measurement precludes a useful com- 
parison with our result; however, they al so calculate the 
mass from the evolutionary m odels of |Burrows et al. 
( |1997[ ) and Chabrier et al. ( 2000 1 . Based on an interpola- 
tion over a grid ot temperature, luminosity, mass, and age 
provi ded by the evolutiona ry models (the Tucson models 
fromlBurrows et al.H1997l and the DUSTY models from 



Chabrier et al. 2000 ) they conclude M tot to be 0.08±0.06, 



consistent with our derived value f r om th e younger age 
range. However, Konopacky et al. (20101 also conclude 



that objects of spectral type late 1V1 through mid L tend 
to have their masses un der-predicted by cur rent evolu- 



tionary models (see also Dupuy et al 



2008). With re- 



fined uncertainties on the mass in combination with the 
age information that can be gathered from NLTT 20346, 
2MASS J0850+1057 will be a useful system for testing 
this result. 

We also calculated the component m asses of NLTT 
20346 using the evolutionary tracks of |Burrows et al.| 
(119971). Using the Lfr Q ; values for an M5 and M6 from 
Keid & Hawleyl (|2005l) combined with the age described 



above, yields component masses of 0.16 ±0.01 M and 
0.13±0.01 M or a total mass for NLTT 20346 of 



0.29± 0.01 M©. This is in agreement with |Delfosse et al 
(2000) where the masses of an M5 and M6 are calcu- 



latecTto be 0.14±0.1 and 0.11±0.1(assuming a 10% un- 
certainty), respectively, using a fourth degree polynomial 
fit to parallax data. 

Combining the mass results for 2MASS J0850+1057 
using the younger age range with those of NLTT 20346 
yields a total mass for this quadruple of 0.36 ±0.02 M©. 
5. DISCUSSION 

We have added NLTT 20346/2MASS J0850±1057 to 
the growing list of VLM (M tot <0.2 M©) multiples widely 
separated from a more massive companion. Table [3] is a 
selected compilation of these widely-separated multiples 
and demonstrates a range in both mass and separation 
of components. Within this list, NLTT 20346/2MASS 
J0850+1057 has a significantly lower binding energy and 
is one of the lowest total mass triple, quadruple, and/or 
quintuple systems known. Figure [8] shows the binding 
energy versus total mass of systems gathered from the lit- 
eraturaHand demonstrates that this new system falls well 
below the binding energy limitation set by known tight 



low mass ( M tnt < 0.2 M^) multiple s (see Close et al 



2007 , 



2003| and |Burga sser et al 
(TOT 



Pinfield et al 



2003). For comparison, 
use the virial theorem to obtain 



the gravitational binding energy of the Pleiades cluster. 
Their derived value (~ 5.3 x 10 38 erg) is smaller than the 
estimated binding energy for the NLTT 20346/2MASS 
J0850±1057 system (and smaller than the lowest binding 
energy systems found to date) , although it is unclear how 
long this association will remain bound. An interesting 
investigation would be to search for stars in the vicinity 
of NLTT 20346/2MASS J0850±1057 as the hierarchical 
nature of this system combined with the common kine- 
matics, potential youth, and wide separation between the 
M dwarf and brown dwarf systems might be indicative 



7 Stellar companions were gathered from th e cat alogs of|Duquen-| 
|noy fc Mayor]jl991| |Fischer fc Marcy||1992| and |Tokovinin| | i99'7r 
and young UCD companion systems from Kraus ct al. 2005 , 2006 
|Konopacky et al |2007||Luhman et al.|2009[ and |Allers et = aTf i |2009| : 
Details on the field UUL) systems were gathered fro m the Very 
Low Mass B inary Archive \http : / / vlmbina ries.org\ see |Burgasser| 
|et al.|2007| and references therein.) 
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of its own moving group. 
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Fig. 8. — A plot of the binding energy vs. total mass. Objects 
marked with filled circles are tight low-mass systems (typically 
Utot < O.2M and p <20 AU). Wide systems (p >100 AU) con- 
taining a UCD companion are marked as five point stars and those 
highlighted in red also contain a VLM binary (hence at least a 
triple). Those marked as squares are systems containing a tight or 
widely separated UCD with an age < 500 Myr. Objects marked 
by open circles come from stellar companion catalogs. The NLTT 
20346/2MASSJ0850+1057 system is the enlarged five-point star. 
The minimum binding ener gy corresponding to tig ht very low mass 
syste ms is labeled (see |Close et a l. 2007 2003| and |Burgasser et ah] 
[20031 - 



Figure [9] shows the M tot vs. sepa ration for th e same 
com panion systems. The studies of Reid et al. (20011 
and |Burgasser et al.| ( |2003[ ) suggested an empirical limit 
for the stability of VLM multiples based on objects that 
were known at the time of the respective works. Nei- 
ther cut-off seems appropriate for NLTT 20346/2MASS 
J0850+1057 or the collection of slightly more massive 
widely sep arated systems plotte d o n Figure [9] (see dis - 
cussions in Fahcrty et al. 2010 and |Dhital et~al.]|2010[ ). 

2"0"l0| ) found a log-normal limitation on the 



Dhital et al. 

separation of binaries in their catalog of 1342 wide (> 
500 AU) low-mass (majority had M tot > 0.3) systems 
and this cut-off does encompass NLTT 20346/2MASS 
J0850AB and all low mass objects on Figure [9j 
The stability of binar y systems will be a fu nction of 

While 



age and total mass (see |Weinberg et al. 
NLTT 20346/2MASS J0850AB is fit by 



1987J). 
re cut-off lim- 



itation of higher mass companions, it is significantly 
different from other known slightly less massive higher 
order multiples so we investigate one p ossible mecha- 
nism for its format ion. Recent work by Sta matellos fe 



Whitworth ( |2009b[ ) has been successful in accounting 
for widely separated VLM binaries using simulations of 
gravitational fragmentation of massive extended disks. 
In their smoothed particle hydrodynamic (SPH) simu- 
lations, a system with Mdisk=^-7Mq, Rd/5x=400AU, 
M star =0.7 M Q is evolved for up to 20 kyr followed by 




10 3 10 2 10' 

Separation UCD to Star (AU) 



Fig. 9. — A plot of the total mass vs. separa tion. Symbols are 
described in Figure \8\ We have over- plotted the lReid etH1|2001[ | 
curve (center) whicn distinguishe d the cut-off lor tne t ormation 
of wide stellar pairs as well as the |Burgasser et al!] (12003} line (far 
right-which is specific for Mtnt <0.2 Ma held systems) and the log- 
normal limitation found in |Dhital et al.| ( |2010[ ) (far left— for systems 
with M tot > 0.3 Mq). 

an N-body dynamical evolution for up to 200 Kyr. Af- 
ter 12 simulations, 96 stars are formed with brown dwarf 
or low-mass secondaries and among those companions, 9 
are tight VLM multipl es. The characteristics of the triple 
systems were listed in Stamatellos fe Whitworth] ( |2009b| ) 
and one was found to have a total secondary mass of 
0.148 Mq, a close brown dwarf-brown dwarf binary sep- 
aration of ~ 1 AU and a wide separation from the central 
star of 7700 AU. This is slightly more massive but highly 
analogous to our proposed system indicating that gravi- 
tational fragmentation could account for the existence of 
NLTT 20346/2MASS J0850+1057. We found no analogs 
created using different formation mechanisms therefore 
comparisons with the predictions from ejection, etc. are 

not possible at this ti me. 

The simulations of |Stamatellos fe Whitworth (2009a I 
showed a significant population ot low mass {M seC ondary 
< 100 Mj) companions at distances out to 10,000 AU. 
This is in agreement with recent studies that have re- 
vealed a growing number of ultracool dwarfs with sepa- 
rations approaching and in some cas es exceeding 10,000 
AU from a companion (s ee Table 1 in|Faherty et al.|2009 
and r eferences there-in; |Dhital et al.||2010|Zhang et al. 
2010 1 . The caveat in this comparison is that the observed 
sa mple covers primarily field age objec ts (1-5 Gyr) while 



the Stamatellos & Whitworth (2009a) simulations stop 



after only 200 kyr ot dynamical evolution. The existence 
of a population of older, widely separated systems sug- 
gests that dynamical interactions into field ages does not 
d isrupt all systems out to 10 000 AU . 



Stamatellos & Whitworth (2009b) find that the nine 
low-mass binaries formed through gravitational fragmen- 
tation that remain bound to their parent star have high 
eccentricities {<cbin> of 0.7 ± 0.2). They postulate 
that the dynamical interactions which form them, and/or 
subsequent dynamical interactions with other stars that 
have condensed out of the disk, account for the increased 
eccentricity of the orbit. The masses of the primary 
stars in that work are not reported although most are 
described as Sun-like. We can compare the eccentrici- 
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TABLE 3 

Details on Triple/Quadruple Systems Containing at Least Two Ultracool Dwarfs 



Name 


SpT 


SpT 


9 


P 


P 


P 


M 


M 


M 


q a 


BE 


Ref 








(") 


(AU) 


(") 


(AU) 


(Mg) 


(•Mo) 


(-M©) 




10 41 Erg 






Primary 


Secondary 


star-UCD 


star-UCD 


UCD-UCD 


UCD-UCD 


Primary 


Secondary 


Total 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


NLTT 20346 


M5+M6 


L7+L6.5 


248 


7700 


0.16 


4.7 


0.290 


0.070 


0.360 


0.24 


0.37 


1A,1B 


G 171-58 


F8 


L4+L4 


218 


9202 


0.33 


10.2 


1.150 


0.095 


1.245 


0.08 


1.7 


2 


HD 221356 


F8 


M8+L3 


452 


11900 


0.57 


15 


1.020 


0.160 


1.180 


0.16 


1.9 


3 


G 124-62 


dM4.5e 


Ll+Ll 


14 


1496 


0.42 


14.3 


0.210 


0.144 


0.354 


0.69 


2.8 


4 


cps Ind 


K5 


T1+T6 


402 


1460 


0.73 


2.6 


0.670 


0.072 


0.742 


0.11 


4.6 


5 


Gl 417 


G0+G0 


L4.5+L6 


90 


2000 


0.07 


1.5 


0.940 


0.143 


1.083 


0.15 


9.4 


6A,6B 


LP 213-67 


M6.5 


M8+L0 


11 


230 


0.12 


2.8 


0.100 


0.176 


0.276 


1.76 


11 


7A,7B 


GJ 1001 


M4 


L4.5+L4.5 


19 


180 


0.09 


1 


0.250 


0.136 


0.386 


0.54 


26 


8A,8B,8C 


G1337 


G8+K1 


L8+L8/T 


•13 


880 


0.53 


10.9 


1.740 


0.110 


1.850 


0.06 


30 


9A,9B 


HD65216 


G5 


M7+L2 


7 


253 


0.17 


5.9 


0.940 


0.167 


1.107 


0.18 


87 


10 


GJ569 


M2.5 


M9.0+M9.0 


5 


50 


0.10 


0.9 


0.350 


0.126 


0.476 


0.36 


123 


11 


Kelu-1 


L0.5/T7.5 


L3p 






0.29 


5.4 


0.125 


0.055 


0.180 


0.44 


178 


12A,12B,12C 


LHS1070 


M5.5 


M8.5+M9.0 


1 


4 


0.45 


3.4 


0.210 


0.138 


0.348 


0.66 


1134 


13A.13B 


HD 130948 


G2 


L4+L4 


0.4 


7 


0.13 


2.4 


1.030 


0.109 


1.139 


0.11 


2178 


14A.14B 



Referenc es. — lA=This Paper lB=B urga sser et al (2 10) 2 = |Faherty et aT] \2010) 3= |Caballero| l|2007) 4=|Seifahrt et al.| (|2005|l 5=|Scholz 
et a lj p003t 6A=| Klrkpatri ck et al.| J200lj) 6B=|Bouv et al. (2003) VA=|Gizis et~al.rpDDb| 7B = Ulosc et al. (200 31 8A=|Golimowski et al.| l2004) 
K irkpat ricFet al. (1999) 8 C=Ma rtin e t al. | |1999) 9A= j Wilso n et al. ffiiUUTt 9B jjjurgass er ct al. (2U05) 10= MugrauereFal. (2007) lT ^Martir i" 
" 2UUU[ 12A= [ Ruiz et al.| jl9y?| 12B=|Liu fc Leggett] pDU5| 12C= |Stump I et ai.| l |2oU8T l 13A= fCemeTret al.||l994| 13B= |Leinert et"al \^2iMl\ 



ct 



a 



'otter et al.| ] |2U02F 14~B= [Dupuy et alTp|Z0Ug| r 



a q represents the mass ratio 

- - pi l iilu. i y 

ties of known widely separated VLM binaries with the 
simulation predictions although the primary masses are 
most likely smaller. Three of the VLM binary systems 
from Table [3j GJ 569B, HD 130948B, and 2MASSJ 
0850+1057, have dyna mical mass measurements and ec- 
centricities measured ( Dupuy et al.| 2008 Konopacky 
et al. 2010[). The mean eccentricity for these three sys 



terns is 0.4 ± 0.2, falling within la of the simulations. 
Further information on the eccentricity of VLM binaries 
widely separated from a more massive component will 
enhance our understanding of dynamical interactions of 
low-mass hierarchical systems. 

6. SUMMARY 

We have presented evidence that NLTT 20346 is 
an M5+M6 binary, and it forms a comoving wide ( 
^7700 AU projected separation) companion system with 
2MASS J0850+1057 at a distance of ~ 31 pc from the 
Sun. NLTT 20346 has moderate chromospheric emission 
and a statistical age based on H a activity of ^5 - 7 Gyr 
for both components. However, it is also an X-ray source 
and has an age of 250-450 Myr based on X-ray activity. 
2MASS J0850+1057 shows strong Li absorption indicat- 
ing an upper age limit between 0.8-1.5 Gyr, favoring the 
younger age for the primary. Assuming an age of 0.25- 
1.5 Gyr for the quadruple, we calculate a total mass for 
2MASS J0850+1057 of 0.07 ±0.02M o and a total mass 
for NLTT 20346 of 0.29±0.01 M Q . The new age will 
be an important parameter for refini ng dynamical mass 
calcu lations of 2MASS J0850+1057 fKonopacky et al. 
20101 so this object can be used as further evidence for 



or against the conclusion that evolutionary models are 
under-estimating the masses for L dwarfs. 

We report a new parallax measurement for 2MASS 
J0850+1057 based on mu lit-epoch ISPI imag ing. Two 
prior measu rements from Dahn et al. (2002) and Vrba 
et al. ( 2004[ ) were discrepant from one another showing a 



2.4cr difference in parallactic angle and a 3.4cr difference 
in proper motion position angle. Our new measurement 
is too uncertain to differentiate between the two previous 
works but a contaminating object located on our recent 
astromctric follow-up images may provide an explanation 
for the differing distance measurements. 

The binding energy for this new quadruple is a mong 
the lowest known for a wide c ompanion system. |Sta-[ 



matellos & Whitworth 
system to IN LIT' 20346 



2009b| presented an analogous 
J0850+1057 formed by 
the gravitational fragmentation of a massive extended 
disk. In that work, they create a significant population 
of wide (p >> 100 AU) systems containing low mass 
(M secon dary < 50 Mj) secondaries analogous to wide 
companion systems now being found. While dynamical 
evolution of their systems ended after just 200 kyr, the 
majority of the wide systems found to date are field aged 
indicating that such systems remain stable for several 
Gyr. 
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